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In this paper we study the superfluid-Mott-insulator phase transition of 

ultracold dilute gas of bosonic atoms in an optical lattice by means of Green 
function method and Bogliubov transformation as well. The superfluid-Mott 
insulator phase transition condition is determined by the energy-band struc- 
ture with an obvious interpretation of the transition mechanism. Moreover 
the superfluid phase is explained explicitly from the energy spectrum derived 
in terms of Bogliubov approach. 
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I. INTRODUCTION 

Ultracold bosons in an optical lattice provide a tunable quantum system with variance 
of the potential depth and lattice constant which can be achieved by adjusting the param- 
eters of the configuration of laser beams. Various quantum phenomena, for instance, Bloch 
oscillations, Wannier-Stack ladders [1] have been investigated in such a system which shares 
spatial periodicity with crystal lattice in solid state physics, however is immune from scat- 
tering of impurities or phonons. The superfluid-Mott-insulator (SMI) phase transition is 
one of the most significant quantum phenomena of condensate bosons in the optical lat- 
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tice. A known analogous system exhibiting the SMI phase transition is liquid helium with 
short-range repulsive interaction in periodic potential [8]. The atomic gas in Bose-Einstein 
condensate (BEC) subjected to the lattice potential which is turned on smoothly can be 
kept in the superfluid phase (SFP) as long as the repulsive interaction between atoms is 
small with respect to the tunnel coupling. With increasing of the potential depth of the 
optical lattice it is getting more and more difficult for bosons to tunnel from one site to 
the other, and finally the system attends an insulator phase above a critical value of the 
potential depth. Considerable attention has been attracted to theoretical researches for un- 
derstanding of the phase transition and determining of the transition condition [2]- [7], in 
which the Bose-Hubbard Model(BHM) is introduced as the starting point of the theoretical 
studies [2] . The phase transition phenomena have been also observed experimentally in BEC 
loaded in a three dimensional optical lattice [2]. Using a strong coupling expansion in terms 
of the hopping term called the decoupling approximation, which is as a matter of fact based 
on the mean-field method D. van Oosten et. al have obtained an analytic phase transition 
condition (see the following Eq.(27)). From an alternative viewpoint the phase transition 
condition can be also determined from the energy spectrum of the system since the excita- 
tion spectrum is neccesarily gapless for the SFP while has a finite gap for the Mott insulator 
phase (MIP). We in the present paper use both the Green function method and Bogliubov 
approach to obtain the explicit excitation energy spectrum and hence to investigate the SMI 
phase transition. 

We begin with the following second-quantized Hamiltonian operator [2] for the system 
of bosonic atoms in the optical lattice, 



where ip (x) and ^ (x) denote the boson field operators which obey the boson commutation 
relation 




(1) 
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(x) , ft (x')j = 5 (x - x') . 



(2) 



Here Vq (x) = X]j=i ^j'o ( x ) sm2 (^r") is the potential of the optical lattice formed by the 



laser light of wavelength A, and hence the lattice constant is d = A/2. Vr(x) denotes 
an external trap potential, and the interparticle interaction is approximated by the short- 
range potential g5(x — x'), where g = A-na s h 2 /m is the coupling constant with a s the s- 
wave scattering length. Expanding the field operator ip (x) in the Wannier basis such that 
-0 (x) = ^ aiW (x — Xj), we obtain the Bose-Hubbard model 



where is the annihilation operator of a particle at the lattice site i, which is assumed as 
being in a state described by the Wannier function w(x. — Xj) of the lowest energy band 
localized on ith site. This implies the assumption that the energy involved in the system is 
small comparing with the excitation energies of the second band. Xj denotes the position of 
the ith local minimum of the optical potential, and nj = a\ai is the number operator. The 
annihilation and creation operators and a\ obey the canonical commutation relations [a* 
,aj ] = 5ij. The parameter J is the hopping matrix element between adjacent sites i, j, and 
is evaluated as 



be of the same value e in the present paper. The interparticle interaction is characterized 
by the parameter 
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For the sake of convenience we rewrite the Hamiltonian(3) in the following form: 





where 

e for i = j 



T 



-J i,j are nearest neighbors 
otherwise 



We see that the first part in Hamiltonian Eq.(6), H , is the same as that of a simple lattice 
under tight-binding approximation (TBA) in solid state physics, so we can rewrite Ty as 

k 

where k is the wave vector in the first Brillouin zone, N s is the total number of the lattice 
sites and e(k) is the energy spectrum of the Hamiltonian H . The inverse transformation is 
written as 

e(k) = N; 1 Tij exp[-ik • (x, - x,)] (8) 
and can be approximated by e(k) w e — ^- exp[—ik ■ (xj — x^)] (TBA energy band) in 

s <ij> 

simple cubic lattice. The explicit energy spectrum is seen to be 

kX 

e{k) = e — Jzcos(— ) (9) 

where z is the number of nearest neighbors of each site. 

The existence of a finite gap in the excitation spectrum is the characteristic of the MIP. 
In section II we attempt to determine the SMI phase transition condition from the energy 
band structure of the ultracold bosonic atoms in optical lattice in terms of Green function 
method. In section III the Bogliubov transformation is used to obtain exact energy spectrum 
with which the superfluid phase is explained explicitly. 

II. GREEN FUNCTION APPROACH 

We begin with the operators, Oj(t), aj(i') , in Heisenberg picture, i.e. a,i(t) = e lHt a,ie~ lHt 
and at(i') = e lHt 'a^e~ tHt ' (in the unit of h = 1). The retarded single-particle Green function 
at zero temperature [9] is defined by 



ai(t); a}(t'))) = -i9(t - {[a^t), $(*)]) (10) 
= -iB(t - t>) [(a^a^t')) - (a](t')a t (t))} 

where 9(t — t') is the step function: 

v ; \ 0, t<? 

The Green function (^(d,i(t); a] (£'))) depends only on the time difference (t — t'). The 
Fourier transformation of the retarded Green function ((z^t)] at(O)')') is seen to be 

G jA (uj) = (laAa) 



- dt ((ai(t); at(0))) exp[i(a; + irj)], i] = +0 



for a real frequency a;. Using Heisenberg equation, we obtain 

uGij{u) = ([a h a]]) + (([a,, H]|aj)) w ( n ) 
which can be evaluated in terms of the commutation relation, 

[a,, H] = ^Tijaj + Uriiai , (12) 
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where () denotes the ground state expectation value. The result is 

uGij(u) = Sij + ^2T im G mj (uj) + UTijiuj), (13) 

m 

where 

Tijfa) = ^Wiai|at)) 
is the higher order Green function which satisfies the following equation 

uFijiu) = ([riiai, at]) + (([mai, #]|a$)) • (14) 

It is easy to find 

[nidi, H] = enidi + ^ T^n^j + ^ T^ajaj - at a*) a* + Unfiidi (15) 



Substituting Eq.(15) into Eq.(14), we see that the obtained equation is not closed, because 
more higher order Green functions appear in the formula of T^ui). We in this stage use 
mean field approximation for the underlined operators in Eq.(15), i.e. the number operator 
Hi is replaced by its expectation value (n,) 

3 3 3 

In the above approximation we have assumed that the average occupation number of Bose 
atoms condensed on ground state in each site of the optical lattice is the same, i.e., (rij) = n . 

Since the obvious symmetry of e(k) = e(—k), we have the equality = Tji according 
to Eq.(8). Utilizing translation symmetry of the Bose system we moreover obtain 

((ofo-)-^)) =0. (17) 

3 

The Eq. (15) is then simplified as 

[riiai, H] « (e + Un )niai + n ^ T^aj. 

3\+i) 

A closed equation for the Green function is seen to be 

LoYi^uj) = 2n 5 i:j + n ^2 T imGij(u) + (e + Un )r i:j (uj) (18) 

m 

from which we find 

= u _ e- Un {25t3 + S T im G m j(u)). (19) 

Substitution of IV,- (cj) in Eq. (19) into (13) yields 

uGij(u) = S i:i + ^2T im G mj (u) H UU ° TT (2Sij + Y] T im G mj (u)) (20) 

uj — e — uno ~f 

which is the equation for the site space Green function Gij{uj). The Green function Gij(u) is 
a function of the position difference (xj — x^) of two sites only for a system with translational 
invariance. The equation (20) for the Green function Gij (uj) can be solved with the Fourier 
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transformation. To this end we express the site space operator Oj in terms of the wave- 
vector-space operator a k 

1 



e ik ' x ' 



-k=Y.' ,kx ^- ( 21 ) 



We can prove that the single particle Green function in the Bloch representation is orthog- 
onal, i.e., 



G kk ,{u) = {\ak\a\, 

= Iv 

N it 



ik-Xi^ik'-Xj //•^•.l-^-t 



e 'e J ( ( a,- a 



where G k {uj) = \\^kWkf/ denotes the orthogonal Green function in Bloch representation. 
The Fourier transformation of the Green function Gij(u) is then seen to be 

^■( w )4s e,k ' (xrxi)Gfc(w) (22) 

Substituting Eq.(22) into Eq.(20) the single particle Green function in the Bloch represen- 
tation is explicitly obtained as 

r(,A = uj-e + Un 

k{ ' (u-e(k))(u-e-Uno)-Uno(e(k)-e) { ' 

Rewrite the solution as the following form 

4(1) A (2) 

where E^ 1 and E k denote the poles of the Green function G k (u), and it is seen that the 
excitation spectrum possesses a band structure such as 

= e, E { k 2) = e{k) + Un (25) 

The lowest band shrinks to a single level of zero band width (see Fig.l). Although the energy 
spectrum Eq.(25) comprises two parts, they may, in a certain case depending on the relative 



values of the interatomic repulsion U and the tunnel coupling J, be merged into the same 
energy band. The energy gap between the two bands is [Fig.l] 

A = E {2) |k=o - = Un - Jz 

When the constant of the interatomic repulsion U is large with respect to the tunnel 
coupling J such that A > 0, a gap exists implying the MIP. With increasing of tunnel 
coupling J the gap width A decreases and finally the two energy bands in the excitation 
spectrum overlap and the gap disappears, indicating the SFP. We then obtain the condition 
of SMI phase transition that 

A = (26) 



, namely, 



zJ n 



which agrees with the result in refs. [7] [8] [10] [11]. 

To se 
(*->0) 



To see the SFP more closely we take the zero wave vector limit of the energy band Ej® 



E {2) ~ e - Jz + Un + ^Jz\ 2 k 2 . 

Under the condition Eq.(27) at which the energy gap between E^> and Ej® disappears, we 
have a gapless Goldstone mode in the excitation spectrum such as 

E exc ~ 1 Jz\ 2 k 2 (28) 

which is different from the result of Bogoliubov theory for the system of weakly interacting 
bosons, in the absence of the periodic potential, where the wave vector dependence of the 
excitation spectrum is linear in the zero wave vector limit so that a non-vanishing velocity 
can exist. Strictly speaking what we obtained here is an ordinary fluid phase. The energy 
spectrum of Eq.(25) determined with the help of Green function method is too simple to 



realize the superfluid phase explicitly. This may be due to the particular procedure of the 
approximation used in the above derivation. It is certainly of interest to study the spectrum 
of bosonic atoms in the BEC trapped in the optical lattice in terms of Bogliubov method to 
see weather or not the system can possess a superfluid phase which we are going to discuss 
in the next section. 

It is worth while to point out that when the interaction between bosons vanishes i.e. 
U — 0, the Green function Eq.(23) reduces to the well known single band solution 

GkH \ u=0 = l — - (29) 

u) — 8 yK) 

for bosons in a periodic potential. 



III. BOGLIUBOV METHOD 

Now we study the energy spectrum of boson atoms in the optical lattice by means of the 
Bogliubov method . Using the relational) the Hamiltonian(6) can be converted to 

H = ^2e{k)a ] k a k + H int 
k 

Hint = ^y" ^ ^k+p,k'+p'0'k' a p' a k a p ■ (30) 
S k,p,k',p' 

Since the number of atoms condensed in the zero-momentum state is much larger than one, 
we have + 1 ~ No 1, where Nq is the total number of condensed atoms. Thus 

we can replace the operator a and a\ by a "c" number V^o- The interacting part of the 
Hamiltonian (30) can be written as (in the order of iV"o) 

H int = ^~ N o + l^f J2 k ^- k + ^-k + 2 «I«fc) 

and the total Hamiltonian is 

UN 2 1 { Un 1 

H = 2ivf + N ° (£ " ZJ) + S fc ( -^(^a-k + a[al k ) + (e(k) + Un )a{a k ^ (31) 
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where ^ fc denotes the sum with exclusion of the term of k = 0. 

The following Bogoliubov transformation is introduced in order to diagonalize the 
Hamiltonian(31) 

h = u k a k + v k a^ k 

(32) 

b k = u k a k + v k a- k . 

We require that the quasi-boson operators b k and b\ satisfy the usual commutation relation, 
[b k , b\] — 1 , which leads to the condition 

u\ -v 2 k = l (33) 

for the coefficients u k and v k and then the Hamiltonian can be written as 

H = E c + Hi + H 2 



where 



E c = ^UN n + N (e - zJ) (34) 



is a constant and 



Hi = J2 k (K + V D + Uno)- 2Un u k Vk)b k bk (35) 

^ = J2 k (^T ( u * + v $ ~ + (W-fc + (36) 

To eliminate the off-diagonal part if 2 we require 

° + ^fc) - (e* + Uno) u k v k = (37) 



where e k = zJ{\ — cos (kd)). Introducing a parameter <f> k such that 

Mfc = cosh^fc v k = sinh0fc 

the conditions (33), (37) lead to the useful relations 
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2u k v k Un 
tanh20 fc = — f = - , rr 
u% + v% s k + U n 

u l + v k = cosh(20 fc ) 



E k 

with which the diagonalized Hamiltonian is obtained as 

H = E C + J2 k E k b{h 
where the energy spectrum E k of quasi-particle is 



E k = ^e k (e k + 2Un ). (38) 

The energy spectrum is different from that of Eq.(25) and is typical for the superfluid. The 
energy gap A g of excitation spectrum is obviously 

E c - N e 



A 



9 N, 



2 

The phase transition condition determined from A g = is 

U 1 



o 

l Un -zJ. (39) 



(40) 



2zJ n 

which shows a factor 2 difference comparing with the condition in Eq.(27). This may be 
caused by the approximation itself. When the energy gap disappears, i.e. A g = 0, the 
dispersion relation of E k {k — > 0) reads 

E k - (zJUn d 2 ) 1/2 k (41) 

indicating explicitly the superfluididy in agreement with the Bogliubov superfluid theory for 
weakly interacting bosons in the absence of the periodic potential. The linear wave vector 
dependence of the excitation spectrum E k ( unlike the ordinary fluid Eq.(28) where E exc 
is proportional to k 2 ) is the characteristic of the sperfluid which gives rise to a persistent 
velocity of superfluid or quasi-particle found as 
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v s =^-§)^={zJUn^y'\ (42) 

For the case of boson atoms with repulsive interaction(a s > 0), the parameters J and U 
are positive, v s is a real number which implies a persistent current. The velocity v s can be 
controlled by the tuning of laser lights which result in the optical lattice. As seen from the 
definitions (4), (5) for J and U, these parameters both depend on the Wannier functions 
which are essentially determined by the potential of optical lattice. Therefore J and U are 
not independently tunable by the adjusting of the laser parameters. In fact when the depth 
of the lattice potential increases, the hopping matrix element J decreases exponentially while 
the matrix- element of the on-site interaction, U , increases. We thereby expect that there 
exists a maximum value of the persistent velocity v s in some particular values of J and U . 

IV. CONCLUSION 

We have studied the Bose-Hubbard model of BEC trapped in a periodic potential in 
terms of Green function method and Bogliubov transformation as well. The condition of 
phase transition between SFP and MIP is determined by the energy band structure of the 
excitation spectrum due to, obviously, the competition between the interatomic repulsion 
and the tunnel coupling. Our result agrees with the condition of SMI phase transition 
obtained in literature. The SFP property of BEC in the optical lattice is explained explicitly 
from energy spectrum derived by means of the Bogliubov approach. It is shown that the 
persistent velocity of the quasi-particle in SFP can be tuned by the adjusting of the laser 
lights which result in the optical lattice. 
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Figure captions 
Fig. 1. Excitation spectrum and energy gap. 
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